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We investigate possible effects of high order nonlinearities on the shapes of infection refugia of 
the Hantavirus epidemic. We replace Fisher-like equations that have been recently used to describe 
Hantavirus spread in mouse populations by generalizations capable of describing Allee effects that 
are a consequence of the high order nonlinearities. We analyze the equations to calculate steady 
state solutions. We study the stability of those solutions under physical conditions and compare 
to the earlier Fisher-like case. We consider spatial modulation of the environment and find that 
unexpected results appear, including a bifurcation that has not been studied before. 



I. INTRODUCTION 



Techniques of nonlinear physics and mathematics are 
finding much application in recent times in biological and 
ecological systems, which, in turn, are enriching non- 
linear science. One example is the spread of epidemics 
0, & @, EL 0] , m particular the Hantavirus [|| that has 
been recently modeled H, H, [H| with the help of the 
Fisher equation [Tl| with internal states representing in- 
fection or its absence in the mouse population. Our pur- 
pose in this paper is to study the consequences of Alice 
mechanisms (to be described below) in the dynamics of 
mice that are the carriers of the Hantavirus. Without 
these effects, the (Fisher-like) equations are 0,11], 
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where we display a 1 — d system for simplicity. Here, 
Ms and Mi are, respectively, densities of susceptible 
and infected mice and M = M$ + Mj is the total 
mouse density, a controls transmission of infection 
on encounter, b and c are rates of birth and death 
respectively and D is the mouse diffusion coefficient. 
The third of the Eqs.{T]) is obtained simply by adding 
the first two and is of the standard Fisher form 0]. 
The first two possess the characteristic features of the 
Hantavirus @ that mice are never born infected and 
that they are unaffected in any other way (for instance, 
they do not die faster) when infected. The logistic 
reaction term in the last of the Eqs.(fT]) is made up of a 
linear growth term (b — c)M and bilinear depletion term 
—M 2 /K. Such a description of population dynamics is 
widespread llj . The steady state homogeneous solutions 
for M are 7]: 0, K (b ~ c). The first of these is unstable 
whereas the second, M = K(b — c) is stable. 



The Fisher analysis of the Hantavirus is based on the 
use of a logistic reaction term, whereas the investigation 
we report in the following assumes a Nagumo term 
which provides an additional zero in the nonlinearity 
relative to the logistic case. The physical content behind 
such a term is the Allee effect, in the presence of which, 
unlike in the logistic case, the zero-M solution is stable. 
If M is small initially, it is attracted to the vanishing 
value; if large, it is attracted to the nonzero value. The 
physical origin of the Allee effect is the possible increase 
of survival fitness as a function of population size for 
low values of the latter. Existence of other members 
of the species may induce individuals to live longer 
whereas low densities may, through loneliness, lead to 
extinction. There is a great deal of evidence for such 
an effect in nature [H, [HI and there have been recent 
reports 14, 15] of theoretical work addressing the effect. 



This paper is set out as follows. The Allee effect is 
described by adding cubic terms to the logistic depen- 
dence and the model is displayed in section 2, both in 
dimensioned and dimensionless forms. The former is im- 
portant to understand the connection of the parameters 
to quantities observed in nature, while the latter facili- 
tates mathematical analysis. We also define two quanti- 
ties, a and %, important to our later development. The 
former is central to the classification of regimes of behav- 
ior and is directly associated with the Allee phenomenon. 
We present a linear stability analysis of the steady state 
solutions in section 3 and argue the importance of the 
threshold a — 0.5 on the basis of a Ginzburg-Landau 
discussion. Restricting attention to the interesting case 
when x exceeds a critical value \c, we carry out numer- 
ical solutions of the model equations in section 4 for the 
simple case of spatially constant a, and in section 5 for 
the richer case of spatially modulated a. Conclusions ap- 
pear in section 6. 
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II. THE MODEL 



Our starting set of equations is, instead of (1), 
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The birth term is quadratic in the mouse density and 
the environment population term is cubic. This is in 
contrast to the Fisher case (1) where these terms are lin- 
ear and quadratic respectively. Needless to say, using a 
quadratic versus linear birth rate has nothing to do with 
sexual versus asexual reproduction. Surely, the mice do 
reproduce sexually. The effective power to be used in 
the equations is the result of a variety of factors influenc- 
ing one another including the probability of encounter of 
mates. The appropriate power of the density can be de- 
termined only phenomenologically, after the fact. Similar 
considerations apply to the quadratic versus cubic envi- 
ronment term. The important distinguishing feature of 
the Allee case is the possibility of extinction that small 
enough populations have to face as a result of the non- 
linearity, and mathematically speaking, the existence of 
two (rather than one) stable fixed points. 

The simplest way to analyze the steady states of the 
system is to first look at the solutions of M from the 
last of Eq. (J2j) . There are three possible homogeneous 
steady state solutions for M: Mq = 0, M mx = (bK + 
\Jb 2 K 2 - IcK)/2 and M mn = (bK - Vb 2 K 2 - AcK)/2. 
By performing the linear stability analysis of these solu- 
tions, we find that Mq and M mx are stable while M mn 
is unstable. However, it is important to note that if 
b 2 < Ac/k, there is only one real solution correspond- 
ing to M = 0. It is convenient to reduce Eq. © to 
dimensionless form by performing the following substitu- 
tions: 



x = XjMlJKD, 

t = M 2 nx T/K, 

m = M/M mx , 

m s = M s /M mx , 

mi = M!/M mx . 



(3) 



We introduce two new quantities that will prove impor- 
tant in the subsequent analysis: 
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the following transformed dimensionless equation set. 
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Using b/M mx = (a + 1)/K and c/M^ x = a/K, we get 



III. LINEAR STABILITY ANALYSIS 

Given that the definition of a in Eq. (JSJ) as a mouse 
density ratio excludes consideration of values of a that 
are not real, the only steady values for m are: 0, a and 
1, and there is only one steady stable state for the total 
population, corresponding to extinction. The solutions 
m = 0, 1 are stable while m = a is unstable. The crucial 
consequence of the Nagumo term (cubic nonlinearity) is 
that the solution m = is stable whereas in the Fisher 
case it is necessarily unstable. Thus, in the presence of 
Allee effects, it is possible that the mouse population 
may vanish completely. Since the system in its steady 
state adopts the values of the stable solutions, we need 
to analyze the mouse density when m — or m = 1. 

The case m — is trivial. When m = 1, we have 
only two valid solutions for (m s ,mi), (1,0) and ((a + 
1)/Xi (x ~ 1 ~ a )/x)- We can observe, by performing a 
linear stability analysis, that when the second solution 
adopts negative values, it is unstable while the first one 
is stable. When x = I + at, there is a transition, the 
first solution becoming unstable and the second becoming 
positive definite and stable. Thus, 1 + a represents the 
critical value of x : 

Xc = 1 + a. 

The results can be summarized as follows: 

(0,0), (1,0) if a> X -l 

{m s ,mi) = <( (0,0), 

((a + l)/x,(X-l-a)/x) if a< X -l 

(7) 

These comments about the fixed points shed light on 
the main issue of interest in this paper, viz. the steady 
state spatial density profiles of mice. Here, we would like 
to stress that in the presence of Allee effects, the infected 
mouse density vanishes for a < b. In their absence on 
the other hand, when the Fisher equation is appropriate 
[7J, the condition involves not only a and b but also c 
and K. As x crosses the value Xo there is a change in 
the stable character of the solutions. The system shows 
an imperfect pitchfork bifurcation: the number of fixed 
points is two for x < Xc and three for X > Xc- The 
bifurcation is imperfect as the system under study is 
not symmetric under reflection, i.e., m s — > — m s and 



3 



rrii — > —mi. By contrast, the Fisher modeling of the 
Hantavirus epidemic shows a transcritical bifurcation [?J . 

Our study of this system shows that a = 0.5 marks 
an interesting threshold that separates two regimes cor- 
responding to larger or smaller values of a. To see this, 
recast Eq.® in the Ginzburg-Landau form, 
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is given by 
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The dynamics is determined by T in that the system 
evolves such that ^jf < 0. The expression for T contains 
two terms. The first term involving the derivative tries to 
minimize T by minimizing the fluctuations in the density. 
The minima of second term, i.e., F(m) are governed by 
the reaction term, f(m) = m(l — m)(m — a) and so we 
have, 

F(m) = to 4 / 4 + am 2 /2 - (1 + a)m 3 /3. 



evolves to a steady state characterized by a homogeneous 
solution, with m = 1 and a nonzero value for the infected 
mouse density, rrii = (% — 1 — a)/Xt m the entire spatial 
domain. In other words, the system evolves towards the 
stable state corresponding to m = 1, and not towards 
the other stable state for m = 0. Such a possibility of 
choice between attractors is a feature associated with the 
Allee effect. It does not arise in the earlier description of 
Hantavirus spread 0] via Fisher-like equations because 
there there is, in that description, only one attractor. If, 
on the contrary, a > 0.5, the total mouse density m van- 
ishes completely in the whole spatial domain. This is, in 
fact, an Allee effect. The observed effect is in contrast 
to the result obtained using Fisher-like equation where 
population can never go to zero because that would cor- 
respond to an unstable state. 

The case a = 0.5 presents a much richer and complex 
behavior of the system. This situation is however, not 
of biological interest because it corresponds to a set of 
parameters of measure zero. Although we have obtained 
some interesting results for that case, we do not display 
them in detail because of their lack of relevance to obser- 
vational matters and comment on them only in passing 
in the Conclusion. 



INHOMOGENEOUS ENVIRONMENT 
(SPATIALLY VARYING a) 



The significance of the threshold value a = 0.5 is that 
F(0) = F(l), i.e., solutions corresponding to both m = 
and m — 1 are equally stable. When a < 0.5, the solution 
m = 1 is relatively more stable than the solution m = 0. 
The situation is reversed for a > 0.5. 



IV. HOMOGENEOUS ENVIRONMENT 
(CONSTANT a) 

Because infected mouse density is always zero irrespec- 
tive of values of a in the subcritical regime \ < Xc we 
focus our attention below only on steady state mouse 
density profiles for the other, more interesting, supercrit- 
ical regime, x > Xc- In this case, there is a possibility of 
getting nonzero mi and the steady state depends on the 
parameter a as well. In order to find the spatial mouse 
pattern, we solve Eqs.© numerically with a given ini- 
tial spatial distribution of susceptible and infected mouse 
density. 

We consider a bounded domain where % = 2 every- 
where ( supercritical case since a < 1) and use reflective 
boundary conditions. We analyze various values of a and 
spatially varying initial densities. Except for a = 0.5, 
we find that the system evolves towards a homogeneous 
steady state, with values corresponding to the steady 
solutions. If we take a < 0.5 and m s (x) — m,(a;) = 
A\ cos(27rwa;)|, a periodically modulated initial condition 
for the initial population, we observe that the system 



Noteworthy effects appear when Allee effects combine 
with a spatially varying a. To study them, we maintain 
the relationship a < x~ 1 throughout the extension of the 
whole domain, so that we have a supercritical solution, 
but assume that the parameters that go into the making 
of a, viz., the birth parameter 6, the death rate c, and 
the environment parameter K (see Eqs. (©), are spa- 
tially varying and consequently introduce a correspond- 
ing modulation into a. 

To understand transparently the physical meaning of a 
postulated modulation of a as we take below, it is helpful 
to consider, only for illustration purposes, the case when 
the quantity Ac/b 2 K is small, and to expand the a ex- 
pression in Eq. ([5]) in powers of that quantity. Then we 
get 



4c 



(10) 



Any inhomogeneity in the environment, represented by 
a spatial variation of c, b or k would be reflected in a 
spatial modulation of a. A larger death rate, a smaller 
environmental parameter or a smaller birth rate would 
cause a decrease in a. 

We have analyzed the effect of such a spatial varia- 
tion in two forms: first by considering a sinusoidal mod- 
ulation with a characteristic wave number and then by 
considering a less regular behavior of a. First we take 
a(x) — 0.5 + S sin(27rwa;), with B < 0.5 so we stay within 
the region where x > Xc and avoid a < 0. Unexpected 
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FIG. 1: Mean values of the total (solid) and infected (dashed) 
population as a function of w, the wavenumber that describes 
the modulation of a(x) 



behavior emerges. The system undergoes a transition as 
the wave number u> of the modulation of a crosses a given 
critical value, oj c w 10. When oj < u> c , the modulation of 
a induces a modulation in the population. This occurs 
around 0.5 for the total population and around x ~2^ a 
for the infected population. But as u> exceeds u) c , despite 
the modulation of a, the total population of mice adopts 
a homogeneous profile with two different values. There 
are two branches, each one corresponding to the steady 
solutions shown in Eq.([7]) for uj > u c . Together with 
the absence of modulation in the total population we ob- 
serve that the infected population survives throughout 
the whole domain. While the total population is homo- 
geneous, both the infected and susceptible population are 
oscillatory in space and take values different from those 
previously calculated. 

Fig Q] displays graphically the results discussed above. 
The mean value of the total and infected population is 
plotted as a function lu, the wavenumber that describes 
the modulation of a(x). A bifurcation is seen at lo c sw 10: 
the system undergoes a sudden change of behavior. The 
total population adopts a homogeneous profile while the 
infected population now survives even in regions where 
a(x) > 0.5. This can be observed in the two plots of Fig. 
[21 where the spatial profile of both the total and infected 
populations are shown for values of ui above and below 
the value ui c . 

Next we provide an example of a more realistic behav- 
ior of the environment that allows us to show that this 
effect survives even when the modulation is not regular, 
provided regions exist with oscillations of wavenumber 
over the critical value. In Fig [3] we plot a situation pre- 
senting an irregular profile for a (also displayed in the 
figure) and the solutions for both branches. 
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FIG. 2: Spatial profile of the total (solid) and infected 
(dashed) mouse population for (a) uj < uj c and (b) uj > uj c 
(only the non-null branch is shown). Space is plotted on the 
x-axis in arbitrary units. 



VI. CONCLUSIONS 

In summary, we have studied what consequences the 
possible presence of an Allee phenomenon might have 
in a population of mice carrying the Hantavirus epi- 
demic. The Allee features have been incorporated though 
a Nagumo, i.e. cubic, term in the nonlinearity. We have 
analytically solved and examined the fixed points of the 
problem without diffusion, and numerically investigated 
the steady state profiles in the presence of diffusion. Our 
findings are that, first, as expected, in the presence of 
the Allee phenomenon, the mouse population can van- 
ish completely. This is not possible when the underlying 
equation is of the type used earlier We also observe 
a dependence from initial conditions that are not present 
in the earlier analysis. These effects stem from the exis- 
tence of two, rather than one, stable solutions. Formally 
stated, under the effect of the Allee phenomenon, the sys- 
tem exhibits an imperfect pitchfork bifurcation instead of 
the transcritical bifurcation observed in the earlier Fisher 
case. The most relevant result is the one obtained by in- 
troducing environmental spatial inhomogeneities. The 
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FIG. 3: Spatial behavior of the total and infected population 
as a (also displayed) varies in space. Each plot corresponds 
to a different branch. Depicted are a (solid line), the total 
mouse population m (dashed line) and the infected mouse 
population rrii (dotted line). Space is plotted on the x-axis in 
arbitrary units. 



observed effect is better observed in the most abstract 
case, when the spatial modulation is sinusoidal. We have 
found that there is a bifurcation in the behavior of the 
system. The nature of this bifurcation is more evident 
when calculating the mean value of the population den- 
sities. We have shown the existence of a critical value 
of the spatial modulation wavenumber uj at which the 
behavior of the systems completely changes, displaying 
bistable behavior that depends on the initial conditions. 
Later, by taking an hypothetic general situation, we have 
shown how this effect operates. 

For the sake of completeness, we describe briefly the 
analysis for the physically (observationally) unimpor- 
tant but mathematically interesting case when a exactly 
equals 0.5. The system behaves in a rather complex 
way. Again, we start with initial conditions of the form 
m s (x) — rrii(x) = A\ cos(27ro;x)|. The interesting feature 
is that now both steady states are equally stable. The 
basin of attraction of each state is such that the ampli- 
tude of the initial modulated condition, plays a relevant 
role. There are three different regimes characterized by 



m 



m 



1,0 
0,8 
0,6 
0,4 
0,2 
0,0 



1,0 
0,8 
0,6 
0,4 
0,2 
0,0 







100 



200 



300 



400 



500 





a 








i 








f- 




v\ ft 

A A 




A A 




A /' 


/ \ 




/ \ 


\\ / ' 




- \ 


b 
















J - 




\J / 


' " ~ " N \ / 




A A 








Y 




- \ 


c 


f \ 




n 




n 




n 




r 




i , 




1 , 1 




i , i 








\ , i 





1,0 

0,8 
0,6 
0,4 
0,2 
0,0 







100 



200 



300 



400 



500 



X 



FIG. 4: The evolution of the initial density profile (a) for 
m(dashed line) and rrn (solid line) towards the steady state 
(c) for a = 0.5 and a = 0.45. Transient patterns are shown 
in (b). Here x-axis is the space and y-axis corresponds to 
density. Units are arbitrary. 



1) the system reaching the homogeneous state m = 1 and 
m i = (x — 1 — a )/x m the whole spatial domain, 2) the 
total mouse density m vanishing completely in the whole 
spatial domain, 3) a steady state being characterized by 
spatial periodic patterns with m oscillating periodically 
between and 1. We also observe a similar kind of os- 
cillation for infected mouse density between rrii = and 
m i = (x ~ 1 ~ a )lx- The two first cases are like those 
encountered already in our analysis above for a > 0.5, 
while the third is exclusive to a = 0.5 As A, the ampli- 
tude of the initial condition, decreases the system goes 
from case 1 to case 2 but with an intermediate regime 
corresponding to case 3. For a = 0.5, we note the exis- 
tence of two critical values of initial conditions at which 
we observe a transition from oscillating periodic struc- 
tures to homogeneous structures. The first critical value 
corresponds to a transition from oscillating structure to 
a homogeneous pattern with m = 0. The second critical 
value corresponds to transition from the oscillating pat- 
tern to the homogeneous one with m — 1. Thus, in the 
the presence of Allee effects, the mouse density depends 
not only on system parameters but also on the initial dis- 
tribution, a dependence that is quite impossible in Fisher 
equation treatments Q. 

Whether rodents of the kind we are describing do or do 
not exhibit these various Allee effects we have described is 
a question for the field biologist to pursue. We hope that 
the interesting consequences that we have predicted the- 
oretically in this paper will stimulate observational work 
in this direction. Work related to the present study (but 
quite different in spirit as well as detail) may be found in 
recent papers by Kenkre and Kuperman [14| on bacteria 
in a Petri dish and more recently in the extensive studies 
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of Clerc et al. [15| who report a number of insights into 
Allee effects on pattern formation. Some of the features 
of the present work that distinguish it from those oth- 
ers is the existence of infected and susceptible subclasses 
in the mouse population, from the mathematical point 
of view the corresponding bifurcations that occur in the 
system, and from the physical point of view the relevance 



of the study to the spread of epidemics and formation of 
refugia. In a future publication we will stress a more 
realistic 2-d depiction of the refugia. 
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